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We study quark thermal recombination as a function of energy density during the evolution of 
a heavy-ion collision in a numerical model that reproduces aspects of QCD phenomenology. We 
show that starting with a set of free quarks (or quarks and antiquarks) the probability to form 
colorless clusters of three quarks differs from that to form colorless clusters of quark-antiquark and 
that the former has a sharp jump at a critical energy density whereas the latter transits smoothly 
from the low to the high energy density domains. We interpret this as a quantitative difference in 
the production of baryons and mesons with energy density. We use this approach to compute the 
proton and pion spectra in a Bjorken scenario that incorporates the evolution of these probabilities 
with energy density, and therefore with proper time. From the spectra, we compute the proton to 
pion ratio and compare to data at the highest RHIC energies. We show that for a standard choice 
of parameters, this ratio reaches one, though the maximum is very sensitive to the initial evolution 
proper time. 
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I. INTRODUCTION 

The features of the proton to pion ratio for Au + Au 
collisions at the highest energy at RHIC [l[ have been in- 
terpreted as the emergence of thermal recombination as 
an important mechanism for hadron production in this 
environment. Recall that in p + p collisions, the proton 
to pion ratio as a function of pt remains basically un- 
changed, never exceeding one, for collision energies rang- 
ing from 19.4 GeV at the Tevatron, 44.6 and 52.8 GeV 
at ISR up to 200 GeV at RHIC 0. This stable behavior 
with collision energy can be understood as a manifesta- 
tion of the fact that in p + p collisions, hadron production 
is dominated by parton fragmentation and the proton to 
pion ratio simply reflects the ratio of the corresponding 
fragmentation functions. In contrast, the proton to pion 
ratio in collisions of heavy systems grows with collision 
energy 0, Q| , reaching and even exceeding one in Au + 
Au collisions at y/s NN = 200 GeV for p t ~ 2 GeV/c. 
While it is possible to understand the different behaviors 
of the above ratios in collisions of heavy systems as par- 
tially arising from the increase of radial flow with collision 
energy 0, [y] , another important ingredient that needs to 
be understood is how these ratios are influenced by the 
relative abundance of protons and pions when these par- 
ticles are produced in the collision of heavy systems at 
the highest energies. 

In its simplest form, the recombination scenario ex- 
plains the formation of low to intermediate pt hadrons 
from the bounding of quarks in a densely populated 
phase space, assigning appropriate degeneracy factors 
for mesons and baryons [7], §]. However, it only ac- 
counts for correlations among quarks in momentum space 
but misses describing the corresponding expected correla- 
tion in coordinate space. An implicit assumption is that 



hadronization happens at a single temperature. However, 
it is known that hadronization is not an instantaneous 
process but rather that it spans a window of tempera- 
tures and densities. For instance lattice calculations 
show that the phase transition from a deconfined state 
of quarks and gluons to a hadron gas is, as a function 
of temperature, not sharp. The question that arises is 
to what extent the probability to recombine quarks into 
mesons and baryons depends on density and temperature 
and whether this probability differs for hadrons with two 
and three constituents, that is to say, whether the relative 
population of baryons and mesons can be attributed not 
only to the degeneracy factors but rather to the dynam- 
ical properties of quark clustering in a varying density 
scenario. 

The detailed answer to the above question is a sub- 
ject that belongs to the realm of confinement phenomena 
and thus to the dynamics of non-perturbative QCD. It 
is however possible to also address it in the context of 
numerical models that reproduce aspects of QCD phe- 
nomenology such as quark clustering at low density and 
color deconfinement at high density. One of such models 
is the so called string-flip model [lOj which has proven to 
be successful in the study of quark/hadron matter as a 
function of density [ll], EH E3] • Other approaches toward 
a dynamical description of recombination, in the context 
of fluctuations in heavy ion collisions, have been recently 
formulated in terms of the qMD model [T3| • 

In this work we compute the meson and baryon ther- 
mal spectra from a dynamical model of quark recombi- 
nation. The work is organized as follows: In Sec. UH we 
introduce the formalism to incorporate the probability 
to form clusters of two and three quarks as a function 
of density and temperature in the description of ther- 
mal particle spectra. In Sec. IIIII we present the basics 
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to implement the numerical simulations in the string-flip 
model. In Sec. |TV] we use this model to compute the 
probability to form mesons and baryons from liberated 
quarks. In Sec. |V]we use these spectra to compute the 
proton to pion ratio and compare the results to experi- 
mental data. We finally present our conclusions and give 
an outlook in Sec. IVII 

II. THERMAL SPECTRA 

In general, the invariant transverse momentum distri- 
bution of a given hadron can be written as an integral 
over the freeze-out, space-time hypersurface £ of the rela- 
tivistically invariant phase space particle density F(x, P), 

dN f P-u(x) , 

E d?p =9 L f dE W F(x ' P) ' (1) 

where P is the hadron's momentum and u{x) is a future 
oriented unit four-vector normal to S and g is the degen- 
eracy factor for the hadron which takes care of the spin 
degree of freedom. 

In the recombination model, the phase space parti- 
cle density is taken as the convolution of the product of 
Wigner functions for each hadron's constituent quark at a 
given temperature and the constituent quark wave func- 
tion inside the hadron. For instance, the meson phase 
space distribution is given by 

F M (x,P) = £ / dz\^ b (z)\ 2 w a ( X ,zP + ) 

x w b (x,(l-z)P+), (2) 

where P + is the light-cone momentum, ^^(z) is the 
meson wave function and a, b represent the quantum 
numbers (color, spin, flavor) of the constituent quark 
and antiquark in the meson, respectively. An analogous 
equation can also be written for baryons. When each 
constituent quark's Wigner function is approximated as 
a Boltzmann distribution and momentum conservation 
is used, the product of Wigner functions is given by a 
Boltzmann-like factor that depends only on the light- 
cone momentum of the hadron [8[. For instance, in the 
case of mesons 

W a (x,zP+)w b (x,(l-z)P+) ~ e ^P+/T e -(l- z) P+/T 



In this approximation, the product of parton distribu- 
tions is independent of the parton momentum fraction 
and the integration of the wave function over z is trivially 
found by normalization. There can be corrections from 
a dependence of each constituent quark Wigner function 
on momentum components that are not additive because 
energy is not conserved in this scenario [HI]. The QCD 
dynamics between quarks inside the hadron is encoded 
in the wave function. 



In order to allow for a more realistic dynamical re- 
combination scenario let us take the above description as 
a guide, modifying the ingredients that account for the 
QCD dynamics of parton recombination. Let us assume 
that the phase space occupation can be factorized into 
the product of a term containing the thermal occupation 
number, including the effects of a possible flow velocity, 
and another term containing the system energy density e 
driven probability V(e) of the coalescence of partons into 
a given hadron. We thus write the analog of Eq. as 

F(x,P) = e- p -^/ T P(e) , (4) 

where v(x) is the flow velocity. 

To compute the probability V(e) let us consider a 
model that is able to provide information about the likeli- 
hood of clustering of constituent quarks to form hadrons 
from an effective quark-quark interaction. We will ex- 
plicitly consider the so called string-flip model which we 
proceed to describe. 

III. STRING FLIP MODEL 

The String Flip Model is formulated incorporating a 
many-body quark potential able to confine quarks within 
color-singlet clusters At low densities, the model de- 
scribes a given system of quarks as isolated hadrons while 
at high densities, this system becomes a free Fermi gas 
of quarks. For our purposes, we consider up and down 
flavors and three colors (anticolors) quantum numbers. 
Our approach is very close to that described in Refs. [ll| 
and [121 ] , where we refer the reader for an extensive dis- 
cussion of the model details. For completeness, here we 
describe the basic ingredients. 

A. Many-body potential 

The many-body potential is defined as the optimal 
clustering of quarks into color-singlet objects, where by 
optimal we mean the configuration that minimizes the 
potential energy. In our approach, the interaction be- 
tween quarks is pair-wise. Therefore, the optimal clus- 
tering is achieved by finding the optimal pairing between 
two given sets of quarks of different color for all possible 
color charges. 

Consider, for example, a set of A quarks some with 
color ci and others with color C2, irrespective of flavor, 
in accordance to the the flavor-blindness nature of QCD. 
We define the optimal pairing between c\ and ci quarks 
as: 

A 

V CxC2 =min^2v[r ici ,P(r iC2 )], (5) 

i=l 

where r, Cl is the spatial coordinate of the i-th c\ quark 
and P(ri C2 ) is the coordinate of the mapped j-th C2 
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quark. The minimization procedure is performed over 
all possible Nl permutations of the C\ quarks and the in- 
teraction between quarks (y) is assumed to be harmonic: 

v ( r ici,fjc 2 ) — 2^( ric i — r jc 2 ) 2 ' (6) 

where k is the spring constant. Through this procedure, 
we can distinguish two types of hadrons: 

i) Meson-like. In this case the pairing is imposed to be 
between color and anticolors and the many-body poten- 
tial of the system made up of mesons is given by: 

V, = V BS + V GG + V RR (7) 

where the individual terms are given by Eq. [5] for the 
corresponding colors. R(R), B(B) and G(G) are the la- 
bels for red, blue and green color (anticolor) respectively. 
Note that this potential can only build pairs. 

\\)Bary on-like. In this case the pairing is imposed to 
be between the different colors in all the possible combi- 
nations. In this manner, the many-body potential is: 

V p = V RB + V BG + V RG (8) 

which can build colorless clusters by linking 3(RBG), 
6(RBGRBG),... etc., quarks. Since the interaction is 
pair-wise, the 3-quark clusters are of the delta (triangu- 
lar) shape. 

According to QCD phenomenology, the formed 
hadrons should interact weakly due to the short-range 
nature of the hadron-hadron interaction. This is partially 
accomplished by the possibility of a quark flipping from 
one cluster to another. At high energy density, asymp- 
totic freedom demands that quarks must interact weakly. 
This behavior is obtained once the average inter-quark 
separation is smaller than the typical confining scale. 

The Hamiltonian for the system built up of N quarks 
of mass mi and momentum pi is given by: 

N 2 

H = '£^ : + V(x 1 ,...,x N ). (9) 

i— 1 1 

In this work, as a first approach, we will study the meson 
and baryon like hadrons independently. Therefore, V = 
V n or V p , depending on the type of hadrons we wish to 
describe. 



B. The variational wave function 

We use a variational Monte Carlo approach to describe 
the evolution of a system of N quarks as a function of 
the particle density. We consider the quarks moving in a 
three-dimensional box whose sides has length a and the 
system described by a variational wave function of the 
form: 

* A (x 1 ,...,x^) = e- A ^ Xl -" x ^$ FG (x 1 ,..,x JV ), (10) 



where A is the single variational parameter, 
V(xi,...,xj\r) is the many-body potential defined in 
Eqs. ([7]) and (JSJ) for mesons and baryons respectively, 
and $re(xi,...,x/() is the Fermi-gas wave function 
given by a product of Slater determinants, one for each 
color-flavor combination of quarks. These are built up 
from single-particle wave functions describing a free 
particle in a box [!§]■ The selection of the variational 
wave function rests in the fact that we are interested 
in the evolution of the system from low to high energy 
densities. The exponential term is responsible of the 
clustering correlations. At low energy density, the 
system is formed by isolated color-singlet hadrons and 
quarks strongly interact inside each cluster; in this 
case, the exponential term of the wave function has a 
big contribution since the average interquark distance 
is of the order of the confining scale. In contrast, at 
high energy density, where asymptotic freedom takes 
place, the interaction between quarks is weak and the 
system looks like a Fermi gas of quarks. In this case, 
the inter-quark separation is much smaller than the 
confining scale and the exponential term effect vanishes. 

As we will show below, the variational parameter 
changes from a fixed value at low energy density (isolated 
clusters) to zero at high energy density (Fermi gas). 

C. Variational Monte Carlo Calculations 

We need to evaluate the expectation value of the 
Hamiltonian operator given by Eq. ^ and minimize it 
with respect to the variational parameter A 

°<*Alfl*>> =0 , (11) 
O A 

In order to do this, we can take advantage of the struc- 
ture of the wave function in Eq. ([TP]) , which allows some 
simplifications in the calculations. After an integration 
by parts, the expectation value of the Hamiltonian oper- 
ator for a given particle density can be evaluated as: 

<H>a=T fg + 2A 2 <VF>a + <^>a, (12) 

where Tfg is the kinetic energy of a free Fermi gas, 
< V >x is the potential energy as defined by Eqs. ([7]) 
and ([8]) and < W >\ is the term that indicates how the 
kinetic energy increases due to clustering correlations. 
This last is expressed as: 

- 1 

<W> X =J2— (x.-y,) 2 , (13) 

— TTli 
i 

where the sum is over all quarks in the system and rep- 
resents the average position of the two quarks connected 
to the i-th quark for a baryon-like hadron or the corre- 
sponding antiquark position for a meson-like hadron. 

A very important property of having written the ex- 
pectation value < H >\ as in Eqs. (fT!?)) and (fTB"]) is that 



4 



2 



1.95 - 

A 

I 1.9 - 

V 

1.85 




1.8 1 1 1 1 1 1 1 1 

0.35 0.45 0.55 0.65 

X 

FIG. 1: (Color on line) Energy per particle from the simula- 
tion (symbols) as a function of the variational parameter at 
low density p = 0.17po for the baryon system. The solid line 
is a fit to the data. 

the two functions to be evaluated (V and W) are local; 
consequently, their expectation values can be calculated 
using Monte Carlo techniques, in particular we use the 
Metropolis method to do the sampling. Finally, this pro- 
cedure is carried out for different values of the particle 
density. 

The variational parameter has definite values for the 
extreme density cases. At very low density it must cor- 
respond to the wave function solution of an isolated 
hadron. For example, the non-relativistic quark model 
for a hadron consisting of 2 and 3 quarks, bound by a 
harmonic potential, predicts, in units where k = m = 1 
that Xjr — > Aott = \/l/2 and \ p — > Ao p = -y/1/3 respec- 
tively; at very high densities the value of A must vanish 
for both cases. In the appendix, we show how to convert 
from the system of units where k = m = 1 to physical 
units. 



IV. PROBABILITIES 

All the results we present here come from simulation 
done with 384 particles, 192 quarks and 192 antiquarks, 
corresponding to having 32 u (u) plus 32 d (d) quarks 
(antiquarks) in the three color charges (anti-charges). 

To determine the variational parameter as a function 
of density we proceed as follows: first we select the value 
of the particle density p in the box, which, for a fixed 
number of particles, means changing the box size. Then 
we compute the energy of the system as a function of the 
variational parameter using the Monte Carlo Method de- 
scribed in the previous section. The minimum of the en- 
ergy determines the optimal variational parameter. We 
repeat the procedure for a set of values of the particle 
densities in the region of interest. In fig. [T]we exemplify 
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FIG. 2: (Color on line) Normalized variational parameter as 
a function of the energy density for the baryon (triangle sym- 
bols) and meson (circle symbols) cases. 



the behavior of the energy per particle from the simula- 
tion (symbols) as a function of the variational parameter 
at low density (0.17po) for a baryon system. The solid 
line is a fit to the data. We can observe that the expected 
values from the solution to the 3-quark system with a har- 
monic potential are properly approached. This is taken 
as an indication of the range where the value of the vari- 
ational parameter should lay for an isolated hadron. 

In order to check whether the variational approach de- 
scribes the expected results in the low and high density 
limits, we display in fig. [2] the energy density dependence 
of the variational parameter A for the baryon, X p , (tri- 
angle symbols) and meson, A w , (circle symbols) cases. 
In both cases, the behavior of the optimal parameter is 
the expected one, as described in detail in the previous 
section. At low densities, A converges to the theoretical 
value (A 07r = 2 -1 / 2 and Ao p = 3 -1 ' 2 ) for baryons and 
mesons respectively. We use these to normalize the data 
such that the values for low energy densities is around 1; 
at high densities the values of both A p and A w approaches 
zero asymptotically. Nevertheless, there are differences 
between the two cases. In the case of baryons, there 
is a deep drop around 0.7 GeV/fm 3 indicating that the 
length scale for quark confinement has increased in the 
medium (thi s unexpected behavior was also reported in 
Refs. [HI, |T3|), to subsequently evolve into the Fermi gas 
domain. For a meson system, the variational parameter 
has a smooth drop. The above indicates a qualitative 
difference in the formation of baryon and meson clusters, 
as a function of the energy density. 

The information contained in the variational parame- 
ter is global, in the sense that it only gives an approx- 
imate idea about the average size of the inter-particle 
distance at a given density, which is not necessarily the 
same for quarks in a single cluster. This is reflected in 
the behavior of the variational parameter X p for the case 



of baryons which, from fig. [2l goes above 1 for energies 
close to where the sudden drop in the parameter hap- 
pens. We interpret this behavior as as a consequence of 
the procedure we employ to produce colorless clusters for 
baryons, which, as opposed to the case to form mesons, 
allows the formation of clusters with a number of quarks 
greater than 3. When including these latter clusters, the 
information on their size is also contained in A. To cor- 
rect for this, we compute the likelihood to find clusters 
of 3 quarks P3. Recall that for 3N quarks in the system, 
the total number of clusters of 3 quarks that can be made 
is equal to N. However this is not always the case as the 
density changes, given that the potential allows the for- 
mation of clusters with a higher number of quarks. P3 is 
defined as the ratio between the number of clusters of 3 
quarks found at a given density, with respect to N. This 
is displayed in fig. [3] 

Therefore, within our approach, we can define the 
probability of forming a baryon as the product of the 
A/Aop parameter times P3, namely 

V p = A/A 0p x P 3 . (14) 

For the case of mesons, since the procedure only takes 
into account the formation of colorless quark-antiquark 
pairs, we simply define the probability of forming a meson 
as the value of the corresponding normalized variational 
parameter, namely 

= A/Ao.. (15) 

The probabilities V p and V v as a function of the energy 
density are displayed in fig. [4] Notice the qualitative 
differences between these probabilities. In the case of 
baryons, the sudden drop found in the behavior of the 
variational parameter is preserved at an energy density 
around e = 0.7 GeV/fm 3 whereas in the case of mesons, 
this probability is smooth, indicating a difference in the 
production of baryons and mesons with energy density. 

These are the probabilities we set originally ourselves 
to look for and that we use in sec El to compute the 
proton to pion ratio. Before proceeding to that analysis, 
it is also instructive to learn more about the behavior of 
other properties of the quark system with energy density. 
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FIG. 3: (Color on line) Percentage of clusters of 3 quarks for 
different energy densities. 
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FIG. 4: (Color on line) Probabilities to form baryons and 
mesons as a function of energy density. 



A. Clusters' size 

To obtain a more precise information about the clus- 
ters' size, we can look at the distribution of particles into 
a cluster, characterized by a mean square radius (MSR). 
To do so, we proceed to calculate the frequency distri- 
bution of clusters with a given value of the MSR. For 
the baryon system, this corresponds to first identify the 
clusters made of 3 quarks -given that clusters with more 
than 3 quarks are allowed to form- and then to compute 
their MSR. For the meson system, as we allowed only 
the formation of clusters with 2 quarks, we compute the 
MSR for all the quark pairs as selected by the optimal 
pairing search. 



The MSR distributions for a set of energy densities 
are shown in figs. [5] and [6] for the baryon and meson sys- 
tems, respectively. They have been normalized such that 
their corresponding area under the curve is one. It can 
be noticed, in both cases, that the distribution becomes 
narrower as the density grows. The physical meaning is 
that as the system gets denser, quarks are brought closer 
to each other. This is precisely the expected behavior 
when the QGP is forming. 

The overlap between the distributions of particle clus- 
ters MSR at high densities and the distribution at low 
density (we take po = 0.2 GeV/fm 3 as reference), indi- 
cates how much the distribution at a given density looks 
like the distribution at the lowest density. This over- 
lap, as a function of density, is shown in fig. [7] for both 
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baryons and mesons. Notice that in the interval 0.35 
< p < 1 GeV/fm 3 , this overlap in the case of baryons 
is larger than the overlap in the case of mesons while in 
the interval 1 < p < 1.4 GeV/fm 3 the behavior is the 
opposite. We interpret this as a further indication of a 
difference in the production of baryons and mesons as the 
phase transition is taking place. 



B. Correlation functions 

In order to further explore the meaning of this re- 
sult, we have computed the two-body correlation func- 
tion, which measures the probability of finding two given 
quarks at a relative distance r. This function is defined 
as [3 
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FIG. 5: (Color on line) Normalized MSR frequency distribu- 
tions for baryons for a set of energy densities. 



where a, (3 denote the collection of all internal quantum 
numbers, such as color and flavor. The importance of 
this observable is that its behavior can tell us if there is 
any correlation between quarks without relaying to the 
pairing information. For example, the two-body corre- 
lation function between quarks of different colors at low 
density, where quarks are in isolated hadrons, will show 
an increase for distances around the hadron size while 
for high densities it won't be able to identify a partic- 
ular distribution and therefore it will become flat. It is 
worth to mention that if we compute the correlation be- 
tween quarks of the same color, the shape will reflect the 
statistical correlation between fermions by a dip at short 
distances [T3 |. Figures [5] and [5] show the correlation func- 
tion between quarks of different color for a set of energy 
densities, for a baryon and meson systems, respectively. 
It is common to normalize the correlation function such 
that the long distance behavior converges to one. How- 
ever, we have chosen the normalization such that the 
uncorrelated long distance behavior converges to zero. 

To quantify the changes in the shape of the correlation 
function we define the function P g , as 
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FIG. 6: (Color on line) Normalized MSR frequency distribu- 
tions for mesons for a set of energy densities. 
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J g(r, p)dr 
Jg{r, p )dr 



(17) 



where we have again chosen the normal nuclear matter 
density po as the reference. We display in fig. [TU] the 
behavior of such function for mesons (circle symbols) and 
baryons (triangle symbols). Both systems exhibit a rapid 
evolution to low values well before the energy density 
reaches the region where the variational parameter drops 
down, although there is a small change for the baryon 
system in that energy (around 0.7 GeV/fm 3 ), and then 
they evolve to become practically negligible around 1.5 
GeV/fm 3 . 



V. PROTON TO PION RATIO 

In order to quantify how the different probabilities to 
produce sets of three quarks (protons) as compared to 
sets of two quarks (pions) affect these particle's yields 
as the energy density changes during hadronization, we 
need to resort to a model for the space-time evolution 
of the collision. For the present purposes, we will omit 
describing the effect of radial flow and take Bjorken's sce- 
nario which incorporates the fact that initially, expansion 
is longitudinal, that is, along the beam direction which 
we take as the z axis. In this 1+1 expansion scenario, the 
relation between the temperature T and the 1+1 proper- 



7 




e (GeV/fm 3 ) 



3 




-0.5 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 

r (fm) 



FIG. 7: (Color on line) Overlap function describing how much 
the clusters (mesons and baryons), resemble the MSR distri- 
bution of the one computed at normal nuclear energy density. 



FIG. 8: (Color on line) Baryon correlation function for a set 
of energy densities. 



time r is given by 



T = Tn [ — 



where r 



Vt 2 



(18) 



Equation (fl~8|) assumes that the 



speed of sound v s changes slowly with temperature. 
A lattice estimate of the speed of sound in quenched 
QCD [13 shows that v 2 increases monotonically from 
about half the ideal gas limit for T > 1.5T C and ap- 
proaches this limit only for T > 4T C , where T c is the 
critical temperature for the phase transition. No reliable 
lattice results exist for the value of the speed of sound 
in the hadronic phase though general arguments indicate 
that the equation of state might become stiffcr below 
T c and eventually softens as the temperature approaches 
zero [3]. For the ease of the argument, here we take v s 
as a constant equal to the ideal gas limit v 2 s = 1/3. 

We also consider that hadronization takes place on hy- 
persurfaces £ characterized by a constant value of r and 
therefore 



where 



dTi = rp dp d(f> dr/, 



1, t + z 

7? = — In , 

' 2 t-z 



(19) 



(20) 



is the spatial rapidity and p, <f> are the polar transverse 
coordinates. Thus, the transverse spectrum for a hadron 
species H is given as the average over the hadronization 
interval of the right hand-side of Eq. (Q]) , namely 



jH 



d 3 P At 



ds 127F F " ( "' P) ' (21) 



where At 



Tf - t . 



To find the relation between the energy density e - 
that the probability V depends upon- and T, we resort 
to lattice simulations. For the case of two flavors, a fair 
representation of the data Q is given by the analytic 
expression 



e/T 4 



1 + tanh 



T-T c 

bT r 



(22) 



with a = 4.82 and b = 0.132. We take T c = 175 MeV. 

For a purely longitudinal expansion, the flow four- 
velocity vector w M and the normal to the freeze-out hy- 
persurfaces of constant r, it^ 1 , coincide and are given by 



= vf i = (cosh r?, 0, 0, sinhr;), 

therefore, the products P-u and P-v appearing in Eq. 
can be written as 



(23) 

m 



p ■ 



P ■ u — m t cosh(?7 — y), 



(24) 



where nit = 
hadron and 



H 



is the transverse mass of the 



V 



2 \E-P; 



(25) 



is the rapidity. 

Considering the situation of central collisions, we can 
assume that there is no dependence of the particle yield 
on the transverse polar coordinates. Integration over 
these variables gives thus the transverse overlap area of 
the colliding nuclei A, namely 



pdpdcj) — * A. 



(26) 



Looking only at the case of central rapidity, y = 0, 
integration over the spatial rapidity variable can also be 
performed straightforward. The integral involved is 



/oo 



u = 2mtKi(x), 



(27) 
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FIG. 9: (Color on line) Meson correlation function for a set 
of energy densities. 




e (GeV/fm 3 ) 

FIG. 10: (Color on line) Normalized correlation function for 
mesons and baryons as a function of the energy density. See 
text for details. 



where x = m t /T and K\ is a Bessel function of the second 
kind and first order. Gathering the above elements, the 
final expression for the hadron's transverse distribution 
is given by 



E 



dN 



H 



2m f A 



d 3 P (2tt) 3 At 



drrKi 



m t 



T(r) 



P[e(r)].(28) 



Figure [TT] shows an example of the the pion and pro- 
ton distributions obtained by means of Eq. (|28p. where 
we have used the values tq = 0.75 fm and r/ = 3.5 
fm and an initial temperature To = 200 MeV. From 
Eq. (fT8|) . this corresponds to a final freeze-out temper- 
ature of ~ 120 MeV. For protons we take a degeneracy 
factor g = 2 whereas for pions g = 1, to account for the 
spin degrees of freedom. Figure [T^] shows the proton to 
pion ratio for three different values of the initial evolu- 
tion proper time r = 0.5, 0.75 and 1 fm and the same 
finial freeze-out proper-time r/ = 3.5 fm, compared to 
data for this ratio for Au + Au collisions at ^Jsnn = 200 
GeV from PHENIX [1]. We notice that the maximum 
height reached by this ratio is sensitive to the choice of 
the initial evolution time. We also notice that the pt 
value for which the maximum is reached is displaced to 
larger values than what the experimental values indicate. 
This result is to be expected since the model assumptions 
leading to Eq. (f2"5| do not include the effects of radial flow 
that, for a common flow velocity, are known to be larger 
for protons than for pions, and which will produce the 
displacement of the ratio toward lower p t values. 



VI. SUMMARY AND CONCLUSIONS 

In conclusion, we have used the string-flip model to in- 
troduce a dynamical quark recombination scenario that 
accounts for the evolution of the probability to form 



a meson or a baryon as a function of the energy den- 
sity during the collision of a heavy-ion system. Wc 
have used the model variational parameter as a mea- 
sure of the probability to form colorless clusters of three 
quarks (baryons) or of quark-antiquark (mesons). We 
have shown that these probabilities differ; whereas the 
probability to form a pion transits smoothly from the 
high to the low energy density domains, the probabil- 
ity to form a baryon changes abruptly at a given critical 
energy density. Within our approach, we attribute this 
difference to the way the energy is distributed during the 
formation of clusters: whereas for mesons the clustering 
happens only for quark-antiquark pairs, for baryons the 
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LU 
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FIG. 11: (Color on line) Proton and pion invariant distribu- 
tions as a function of transverse momentum for to = 0.75 
fm, Tf = 3.5 fm and To = 200 MeV, corresponding to a final 
freeze-out temperature of ~ 120 MeV. 
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FIG. 12: (Color on line) Proton to pion ratio as a function of 
transverse momentum for three different values of the initial 
evolution proper-time to = 0.5, 0.75 and 1 fm and the same 
finial freeze-out proper-time r/ = 3.5 fm, compared to data 
for this ratio for Au + Au collisions at t^/snn = 200 GeV 
from PHENIX. The height of this ratio is very sensitive to 
the choice of the initial evolution time. 

energy can be minimized by also forming sets of three, 
six, etc., quarks in (colorless) clusters. These produces 
competing minima in the energy that do not reach each 
other smoothly. We interpret this behavior as a signal for 
a qualitative difference in the probability to form mesons 
and a baryons during the collision evolution. 

Notice that this approach is quantitatively different 
to the conventional coalescence scenario [7j where the 
probability to form a bound state of quarks is simply 
given by the overlap of the quark distribution functions 
with the Wigner function of the formed hadron and thus 
ignores the underlying evolving dynamics of the collision 
process. 

We have incorporated these different probabilities to 
compute the proton and pion spectra in a thermal model 
for a Bjorken-like scenario. We use these spectra to com- 
pute the proton to pion ratio as a function of transverse 
momentum and compare to experimental data at the 



highest RHIC energies. We argue that the ratio com- 
puted from the model is able to reach a height similar 
to the one shown by data, although the maximum is dis- 
placed to larger p t values. This could be understood 
by recalling that the model does not include the effects 
of radial flow which is known to be stronger for protons 
(higher mass particles) than pions. The inclusion of these 
effects is the subject of current research that will be re- 
ported elsewhere. 

Appendix: Physical units 

The simulation was performed taking m = k = 1 . Here 
we show the conversion to physical units. 

Baryons: To fix the the energy unit we first notice that 
in a 3-body system the energy per particle, including its 
mass, is given by (with m = k = 1): 

| = V3 + 1. (29) 

If we identify the state as the proton of mass M p = 938 
MeV, then the correspondence is 

V3 + 1 -> 312.7 MeV. (30) 

To fix the length unit we use t he mean square radius, 
which for a 3-body system is: V < r 2 > — (3) 1 / 4 . The 
experimental value for the proton is 

V< r 2 > = 0.880 ±0.015 fm. (31) 

Then the correspondence is: (3) 1 / 4 — > 0.88 fm. 

Mesons: In a similar fashion we obtain for mesons (tak- 
ing the pion as the representative 2-body particle): En- 
ergy: ^ + 1 -► 70 MeV, length: 2 1 / 4 -> 0.764 fm. 
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